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Abstract 



We develop a new symmetrized version of the Projector Quantum Monte Carlo 
(SPQMC) method which preserves the symmetries of the system by simultane- 
ously sampling all symmetry related Ising configurations at each MC step and 
use it to study the effect of electron correlations on a single Ceo molecule and 
its structural motifs, within the Hubbard model. Thi SPQMC method allows 
more accurate estimates of correlation functions as seen from calculations on 
small systems. The method applied to some molecular fragments of Ceo, viz, 
pyracylene, fiuoranthene and corannulene, gives a good chemical description 
of these systems. Analysis of the ground state bond orders allows us to visu- 
alize pyracylene as a naphthalene moiety with weakly bridged ethylenic units, 
fiuoranthene as weakly bridged naphthalene and benzene units, and corannu- 
lene as radialene with ethylenic bridges. We study the ground state properties 
of a single fuUerene molecule, with and without bond alternation. The bond 
orders in the ground state for the two types of nearest neighbour bonds are un- 
equal, even for uniform Ceo. The spin and charge correlations give a consistent 
picture of the interacting ground state in all these systems. We construct pro- 
jections of the bond orders and spin-spin correlation functions on the space of 
irreducible representations of the icosahedral point group. These projections, 
analogous to the structure factors of translationally invariant systems, give the 
amplitudes for distortions which transform as the irreducible representations 
of the point group. The amplitude for the Hg distortion is the largest while 
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the spin structure has large weights in the and Gu representations. 
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1 Introduction 



The fascinating and esoteric phenomena exhibited by systems containing the 
fullerene, CeoiH, 3, have attracted the attention of experimentahsts and the- 
orists in both physics and chemistry. Pure Cqq in the soild state exhibits 
interesting plastic crystal transformations!^. Alkali doped compounds of Ceo 
are superconducting, with i^fei.sCsi.sCgo having the highest known supercon- 
ducting transition temperature for organic systems (T^ = 32K)[0. The donor- 
acceptor compound, Cqq—TDAE (TDAE = tetrakis dimethyl amino ethylene) 
is one of the few known purely organic ferromagnets, with a Curie temperature 
of 16.5K, the highest known for an organic system]^, p. 

Despite intense efforts at understanding the diverse electronic phenomena 
in Cqo systems, a reliable correlated electronic structure for the ground state of 
even an isolated Cqq molecule still remains elusive. A perturbation theoretic 
(PT) calculation in the framework of the Hubbard and modified Hubbard 
models §] focuses on the energies of the low-lying states of doped Ceo- The 
main concern of this study is the sign of the low-spin high-spin energy gap 
as a function of the model parameters. This problem has also been studied 
by the Variational Monte Carlo (VMC) method using a Gutzwiller type 
variational wavef unction [|lO| j . The large U limit of the Hubbard model on a 
single fullerene molecule leads to a frustrated, antiferromagnetically interacting 
s = 1/2 Heisenberg spin model. The classical ground state (large s limit) has 



been determined[ll|. Using a spin- wave analysis around the classical ground 
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state, the possibilty of a magnetic instability in the ground state at a critical 
value Uc of the Hubbard parameter has been conjectured|T2[ . 



In recent years, there have been some reports of calculations of the ground 
state of Cgo using conventional quantum chemical techniques. One of these 
employs an STO — 3G* basis in the Hartree-Fock procedure and follows 
this up with a limited configuration interaction (CI) procedure ||13||. Another 
calculation employs the Modified Neglect of Differential Overlap (MNDO) 
procedure while yet another the Linearized Muffin Tin Orbital (LMTO) 



method |]15|. Even at the level of a single fuUerene molecule, conventional 
quantum chemical techniques for obtaining the correlated electronic structure 
are inadequate. Besides being basis set dependent, these approximations are 
known to be uncontrolled due to arbitrary truncation of the CI procedure. 
Thus, such calculations fail to isolate the role of the active p,r-type orbital in 
determining the properties of fuUerene. Given these difficulties encountered 
in quantum chemical approaches, what is essential to understanding the prop- 
erties of Ceo is a thorough study of well defined model hamiltonians. Such 
studies allow us to follow the evolution of the electronic states as a function 
of model parameters, thereby providing valuable insights. 

The solution of model hamiltonians can also be undertaken within the 
framework of quantum chemical techniques, but even here, the astronomically 
large dimensionality of the full many body space corresponding to the half 
filled (or nearly half filled) p-,r-like orbitals cannot be dealt with reasonably. 
However, exact calculations are possible on chemically feasible fragments of 



Ceo and these studies would provide the necessary cahbration for any novel 
approximate technique p!6|, [1^. 

Amongst recent many body techniques for solving interacting model hamil- 
tonians, Quantum Monte Carlo (QMC) methods and the density matrix renor- 
malization group (DMRG) method llTsj have come to the fore. The DMRG 
method is accurate mainly for one- dimensional lattices. There are many vari- 
ations of the QMC method. While the VMC method [|19|] provides an estimate 
of expectation values of the Hubbard Hamiltonian in the Gutzwiller trial func- 
tion, the temperature QMC methods||2^ extrapolate ground state properties 
from finite temperature properties. A Quantum Monte Carlo procedure that 
directly estimates the ground state properties of Hubbard type models is the 
Projector Quantum Monte Carlo (PQMC) method|21|. In this method the 
ground state is projected from a non-interacting trial wavefunction such as 
a Slater Determinant in the Huckel Molecular Orbital (MO) basis. Also, the 
PQMC method does not suffer from either the limitations of the structure 
imposed by a trial wavefunction, as in the VMC method, or the necessity 
of extrapolating the ground state from finite temperature results, as in the 
temperature Monte Carlo methods (which suffer from numerical difficulties at 



very low temperatures) p2 



In this article, we report results of the first PQMC studies of the ground 
state of an isolated fuUerene molecule. We have greatly enhanced the accuracy 
of the Monte Carlo sampling by exploiting the symmetries of the system. From 
comparison of exact results with symmetrized PQMC results for small clusters, 
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we find that incorporating symmetries leads to a marked improvement in the 
correlation functions, thus providing us a more reliable picture of the correlated 
ground state of Cqq. 

This article is organized as follows. In the next section we describe the 
Projector Quantum Monte Carlo method for the Hubbard model. In section 3, 
we introduce the symmetrized PQMC (SPQMC) method and compare results 
for small systems with results obtained from exact calculations and the VMC 
method. Section 4 pertains to our results on the ground state properties of 
Ceo and its molecular fragments and discussion of the data. Section 5 is a 
summary of the work presented in this paper. 

2 The PQMC Method 

The Monte Carlo algorithm to project out the ground state from a trial 



wavefunction was first introduced by Sorella et a/|^. While their method 
used a fast Fourier transform (FFT) technique for propagating the single- 
particle states and Langevin dynamics to sample the distribution, Imada and 



Hatsugai||23|| introduced an algorithm that runs parallel to the temperature 
QMC, which we employ in our calculations. 

The single band Hubbard Hamiltonian, H may be written as a sum of a 
non-interacting part, Hq and an interaction part. Hi, given by[p4 



H = Ho + Hi (1) 
^0= E T.t^J{cUJa + h.c.) (2) 



<ij> o- 



Hi^U^ fii^fiii (3) 

i 

where (qo-) is the creation(annihilation) operator for an electron with spin 
a in the Wannier orbital at the i*'* site, the corresponding number 

operators and the summation is over bonded atom pairs. 

The ground state, j-^o > of a Hamiltonian, H can be projected from a trial 
wavefunction |0 > by the ansatz 

1^0 >= hm , ''^.^ (4) 

^< (j)\e-^l^^\(f) > 

provided \4> > has a nonzero projection on to \ipo >. The trial wavefunction 
10 > is usually formed from the MOs generated by diagonahzing the non- 
interacting problem. The trial state with M-fermions, in second-quantization 
is represented as 

M N 

\<l>'r>- n(E(*<^U4)|o> (5) 

m=l i=l 

where is an N x M matrix of the MO coefficients whose row index, i, 
corresponds to the site number and column index, m, corresponds to the m*'* 
MO occupied by an electron of spin a, in the trial wavefunction. 

To carry out the projection for an interacting Hamiltonian, it is necessary 
to Trotter decompose the density operator exp{—PH) as 

e-f^^ = (e-^^^)^ (6) 

where (3 — L x At. The error introduced in this decomposition due to non- 
commutativity of and Hi is of order (Ar)'^. The matrix representation, bo, 
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of exp{—AT/2Ho) can be easily obtained in the basis of Wannier functions as 

bo = exp[-K] (8) 

^ _ \ —^tij, for i,j bonded sites; /gv 
1 0, otherwise. 

However, the interaction part cannot be represented as a matrix in the one- 
particle basis. This makes it necessary to obtain a one-particle like description 
of expi^ArHi), which is indeed possible via a discrete Hubbard-Stratanovich 
(H-S) transformation at every time slice. At the Z*'^ time slice, the H-S trans- 
formation gives 

g-AT^i(0 ^ J2 ■■■ J2 e~^''^iT{'.S()g-Ar^u(i,sO ^-^Q^ 

Sll=±l S2(=±l Sjvi=±l 

N -, TT 

^-ArHMi,sO = Uio^M>^^iini^ - Ar^n,T)) (H) 
i=i ^ ^ 

^-ArH,,(i,s,) ^ f^(±^^p^_xs.^n,^ - Ar^n^O) (12) 



.=1^2 



where 



A = 2arctanh^tanh(ATU/A), (13) 

Sii is the Hubbard-Stratanovich field at the i^^ site in the time slice and Si is 
an A^"- vector whose i*'* component corresponds to the Ising variable Su. Now it 
becomes possible to obtain a matrix representation, hia{l, Si), oiexp{—ATHi^{l, si)) 
in the basis of Wannier functions, 

1 At 
(bi|(/, si))ij = dij-exp[Xsu - — ] (14) 
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(bu(/, si))ij = dij-exp[-Xsii - —]. (15) 
This gives us the matrix representation, Bo.(Z,s/), for exp{—ATHa) at the Z*'* 



time shce, 



Ba{l,si) = bobi^(/,s/)bo 



(16) 



The Trotter decomposition followed by the H-S transformation thus results 
in the replacement of the interacting Hamiltonian by a sum of 2^^^ nonin- 
teracting Hamiltonians corresponding to each configuration of the {d + 1)- 
dimensional N x L Ising lattice. The properties of the interacting Hamilto- 
nian are simple averages taken over the corresponding properties of all the 
2NxL noninteracting Hamiltonians. The matrix representation of the projec- 
tion operator for the interacting Hamiltonian can now be written as a sum 
of individual projections, P({s}), each corresponding to the non- interacting 
Hamiltonian for the Ising configuration {s}. The matrix P({s}) is given by. 



This allows the representation of the expectation value, w, of the density 



P({4)=PT({4)Pi(M) 



(17) 



Pa{{s}) = B,(l, Si)B,(2, S2) . . . B,(L, sl) 



(18) 



operator exp{—l3H), in the trial state |0 >, namely. 



w = 



(19) 



which takes the form 



w ^ {<P\^P{{s}M 



(20) 
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within Trotter error. It is convenient to define weights for each Ising configu- 
ration {s} as 

W{{s})^W^{{s})W^{{s}) (21) 

W„{{s}) = det(^lB,{l, Si)B,(2, S2) . . . B,(L, 8^)$,) (22) 

w can now be conveniently expressed as a sum over weights of all Ising con- 
figurations, 

^ = E^({4)- (23) 

For large N and L, it is computationally prohibitive to obtain properties as 
exact averages over all Ising configurations. A simple average based on crude 
sampling would lead to large errors in the estimates. However, by resorting 
to importance sampling of the Ising configurations, it is possible to obtain 
estimates with acceptable errors. A straightforward importance sampling of 
the Ising configurations assumes that an Ising configuration {s} appears with 
probability 

^ ™ (24) 

However, computing w is equivalent to obtaining the ground state projection 
for the given This vicious circle is broken by employing an ergodic Markov 
chain. 

The state space of the Markov chain is taken to be the space of Ising 
configurations and the initial Ising configuration is chosen assuming a uniform 
probability distribution in the configuration space. The one-step transition 
probabilities between any two Ising configurations {s} and {s'} is so chosen 
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that the ratio of the hmiting probabihties of the configurations tends to 
The one-step transition probabihty 

r = ^^^ (26) 
Wi{s}) ^''^ 

gives the correct hmiting ratios for the probabihties of Ising configurations. 
The Monte Carlo simulation now follows the usual procedure. The new Ising 
configurations can be generated either by the cluster spin-fiip or the single spin- 
fiip mechanisms. For convenience we employ sequential spin fiips exhausting 
all sites in a given time slice before moving over to the next time slice. 

In the Monte Carlo procedure, we estimate the equal time Green function 
G„{1, {s}), which we regard as a matrix in the Wannier basis, for spin a, at 
the time slice, for every configuration {s} which is sampled. The (j, fc)*^ 
matrix element of Gcr{l, {s}) is given by 

where the the projected states {L„{1, {s})\ and |-Ro-(/, {s})) are obtained anal- 
ogous to 10 > in eqn. (5) using the coefficient matrices L and R, 

Kil, {s}) = $JB^(1, Si)B,(2, S2) . . . B^(/ - 1, Si„i)bo (28) 

R,(/, {s}) = bi,(/, sOboB,(/ + 1, Si+i)B,(/ + 2, 8^+2) • • • B,(L, 8^)$,. (29) 

In practice, this involves obtaining the inverse of a matrix. For the single spin 
fiip mechanism, efficient algorithms exist to obtain the ratio, r, of the deter- 
minants and to update the Green function if the spin fiip is accepted[^, p5| . 
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However, the updated inverse degrades fairly rapidly and the Green function 
will have to be recomputed from scratch every few hundred spin flips. Prop- 
erties of the Hamiltonian can be evaluated from Monte Carlo estimates of the 
equal time Green function ((G)) 

((^^))^^-j;^T.iGAi,{s}h. (30) 

l,{s} 

where Ng is the number of MC sweeps over the Ising lattice, after equilibration. 
We employ Wick's theorem to compute two particle properties from the single- 
particle Green function, {Ga{l, {s}))ij, at every time slice, /. 

3 Symmetrized PQMC Method 

A serious drawback of Monte Carlo simulations carried out as described in 
the previous section is that the Green function and correlation functions do 
not reflect the symmetries of the system. This is because the noninteracting 
Hamiltonian corresponding to an Ising configuration does not have the full 
symmetry of the system. In the MC sampling each Ising configuration is dealt 
with independently. Thus, when a particular Ising configuration is visited in 
the course of sampling, it is not guaranteed that all the other Ising configura- 
tions related by symmetry are also visited. For example, while we may sample 
the configuration in fig. la, it is not certain that the configuration in fig. lb, 
which is related by C4 symmetry is also sampled during the simulation. For 
this reason, the MC sampling breaks the symmetry of the Hamiltonian. The 
full symmetry of the Hamiltonian can be preserved if we ensure that for every 

13 



Ising configuration that is sampled, all its symmetry related partners are also 
tacitly included in the sampling. 

The symmetry of the Hamiltonian guarantees that the one-step transition 
probability between Ising configurations {s}, {s'} is the same as that between 
R{s} and R{s'}, where is a symmetry element of the Schrodinger group of 
the Hamiltonian, 

P{s}-^{s'} = Pr{s}^r{s'} (31) 

Therefore, if an Ising configuration {s'} is accepted (rejected) from an ini- 
tial configuration {s}, the same result is expected from all symmetry related 
pairs of configurations R{s'} and R{s}. This feature can be incorporated by 
constructing a symmetrized Green function as follows: 

Gri{^})^lT.G^AR{s}) (32) 

R 

and R runs over all the h symmetry elements of the group. 

It appears from the equation that we need to update the Green functions 
G{R{s}) for every symmetry operation R, which could be enormously compute 
intensive. However, the Green functions G{R{s}) and G{{s}) are related as 

G,,{{s'}) ^ G,A{s}) (33) 

R:{s}^ {.'}; (34) 
R-':t^t' ; R-' : j ^ j' (35) 
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Thus, from the Green function of a single Ising configuration, we can generate 
the Green functions of aU Ising configurations related by the symmetry group 
of the system. 

In highly symmetric systems such as Ceo, the number of symmetry ele- 
ments, h, of the Schrodinger group can be very large. For instance, the number 
of symmetry elements in the icosahedral group is 120 and the number of pairs 
of sites is 0{N'^), which, for Ceo is 3600. In a symmetrized sampling procedure, 
for Ceo, it appears that wc should average all the 3600 elements of the single 
particle Green function over the 120 symmetry related Ising configurations. 
This procedure turns out to be prohibitive, computationally. However, we 
recognise that there are only 24 unique pairs of carbon atoms on the fuUerene 
molecule. It is sufficient to average these 24 pairs of matrix elements of the 
Green function over all the symmetry related Ising configurations, to obtain 
the complete symmetrized Green function. This reduces the cpu requirement, 
in the case of Ceo by a factor of 150, thus rendering the SPQMC procedure 
computationally feasible for large systems. 

The PQMC method in its current implementation is not variational. A 
variational implementation of the PQMC procedure can, in principle, be acheived 
by projecting the states {L^{1, {s})\ and \Ra{l, {«})) from identical Ising con- 
figurations. The Green function is sampled only at the middle time slice. This 
ensures that the estimated quantity is an expectation value of a projected trial 
state. However, this version of PQMC would have larger errors for the same 
13 due to halving of the Ising lattice. 
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We have tested the SPQMC procedure for the simple case of a Hubbard 
ring of six sites since this is easily amenable to exact studies. In table 1, we 
compare the SPQMC ground state energies for different values of the Hubbard 
parameter, U/t, with the results of exact and VMC calculations. It is not 
surprising that in terms of energy, the VMC results are better than the PQMC 
results for very small 0.5i), since the Gutzwiller trial function is the true 
ground state lor U/t — 0.0. However, for larger values olU/t, the percentage 
errors in PQMC estimates of the ground state energy are smaller than in the 
VMC estimates. The most striking difference, however, is seen in two-particle 
correlation functions like the spin-spin and charge-charge correlations. The 
results presented in tables 2a and 2b show clearly that the SPQMC results 
are accurate even for U — 4.0t. Thus, the SPQMC method allows a better 
understanding of the correlated ground state of Cqq. 

Symmetrized sampling can also be used to estimate properties of excited 
states. This is possible since the projection operator projects the lowest energy 
state in a chosen symmetry subspace of the Hamiltonian. We could start with 
an initial state that transforms as the basis of an irreducible representation of 
the desired space, usually an appropriate linear combination of Slater Determi- 
nants obtained from group theoretic projection operators. In the symmetrized 
sampling procedure, the symmetry of the Hamiltonian is explicitly preserved. 
This prevents any admixture of states from different symmetry spaces thereby 
leading to the lowest energy state of the subspace. When the starting state 
is a linear combination of Slater Determinants, estimates of properties require 
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modifications to the PQMC procedure as will be discussed elsewhere p6[ . 

For simple cases with non-degenerate LUMOs, such as Hubbard chains the 
properties of the ground state of the doped system can be obtained from a 
single determinantal PQMC with symmetrized sampling, as verified for small 



systems|2g]. The doped states of molecules with degenerate LUMOs are rather 
difficult to handle. To estimate the ground state properties of such systems, 
with open shell electronic configurations, it is necessary to carry out a multi- 
determinantal, symmetrized PQMC simulation. This is true for even as simple 
a system as doped benzene. There is considerable interest in the electronic 
structure of doped Ceo species. Unfortunately, the threefold degeneracy of the 
LUMOs of Cgo makes a single determinantal calculation feasible only for Cqq 
and the high spin state of Cqq . We are currently implementing a multideter- 
minantal, symmetrized PQMC method for studying doped Ceo species. 

4 Results on the Ground State of Ceo and its 
Fragments 

In this section we present results on the ground state of Ceo and its fragments, 
obtained via SPQMC calculations. We compare these with VMC results and 
also with exact results, wherever possible |T^, In Section (4.1) we discuss 



the nature of the correlated ground state in the molecules pyracylene, fluoran- 
thene and corannulene (fig. 2). These molecules have been chosen because they 
represent the characteristic chemical motifs of Ceo and contain one or more 
five-memebered rings essential for displaying the minimal topological features 
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of Ceo- The trends observed in these systems are expected to be indicative of 
the behaviour of the full molecule. In Section 4.2 we present SPQMC results 
on both uniform and bond-alternated Cqq for various values oiU /t. These re- 
sults include energies, bond orders, charge and spin correlations. We compare 
these results with the VMC results. 

4.1 Fragments of Ceo 

In table 3 we present the ground state energies of the fragments, pyracy- 
lene, fiuoranthene and corannulene for different values oiU/t. The pyracylene 
molecule is the smallest characteristic fragment of Cgo that we have studied. 
It is also amenable to exact calculations. The ground state energy obtained 
from the SPQMC calculations compares favourably with that from the exact 
calculation for U /t = 1.0. We also compare our SPQMC results with the VMC 
results for all values of f//t, for all fragments. The agreement is very good with 
VMC for small values of U / 1 and tends to deteriorate for larger values of the 
Hubbard parameter. The SPQMC results are more reliable than the corre- 
sponding VMC results, for similar sample sizes. For, the SPQMC estimates 
are exact within Trotter error (which in our case is 0(10^^)). In contrast, 
the VMC estimates are also limited by the choice of the trial wavefunction 
and it has been shown that the Gutzwiller function, with a single variational 
parameter, is a poor trial function at large U/t p8|| . 

The bond orders in pyracylene, fiuoranthene and corannulene decrease with 
increasing U/t, which is consistent with increased covalency of the ground 



18 



state. A detailed analysis of the bond order data allows one to visualize these 
molecules in terms of interconnected chemical motifs. We can describe pyra- 
cylene as a naphthalene moiety weakly coupled to two ethylene moities. Fluo- 
ranthene can be visualized as a naphthalene motif which is weakly linked to a 
benzene unit. And corannulene can be viewed as ethylene bridges connecting 
the spokes of the radialene molecule. 

The spin correlations of the nearest neighbours in these molecules lend 
support to the picture emerging from the bond order data. The strong bonds 
also have large antiferromagnetic (a/m) correlations. Thus, the strongest bond 
(1-2 bond) in pyracylene also has the largest afm spin-spin correlation. In flu- 
oranthene, the bridging bond (4-5 bond) has a spin-spin correlation which is 
less than a third of that of the strongest bond (7-8 bond). The bonds along 
the perimeter of corannulene fall into two distinct classes, a trend expected 
from the bond order variations. What is more interesting are the bonds of the 
radialene motif (fig. 3), all of which have reasonably large afm correlations. 
However the spokes have larger correlations than the bonds of the pentagon. 
Thus we see that the pentagon in corannulene is already in a chemical envi- 
ronment similar to that of Ceo- Frustration due to the five membered ring is 
thus manifest in a qualitatively different way, in this molecule. 

The charge correlations of the bonds in pyracylene, fluoranthene and coran- 
nulene increase with increasing U/t, while the on-site (diagonal) charge cor- 
relations decrease with increasing U/t, a feature expected of the Hubbard 
model. In the large U/t limit, the intersite charge correlations approach unity. 
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The magnitude of the charge correlations can be related to the total oscillator 
strength of the molecule for dipole transitions This evidently decreases 
with increase in the strength of electron correlations. The diagonal charge cor- 
relations in the noninteracting limit have a value of 1.5 while in the strongly 
interacting case, this valuse is 1.0. However, depending upon the topology of 
the bond connectivities, the interactions lead to values of the diagonal charge 
correlations in the range 0-2. The picture as obtained from charge correla- 
tions computed using the SPQMC method corresponds to a more uniform 
distribution of charges than that obtained from the VMC calculations. 

4.2 SPQMC Ground State Properties of Ceo 

We have computed the ground state properties of an isolated fuUerene molecule 
using the SPQMC method. We have studied uniform fuUerene (equal hexagon- 
hexagon (h-h) and hexagon-pentagon (h-p) transfer integrals) as well as bond 
alternated fuUerene (h-h transfer integral larger in magnitude than the h-p 
transfer integral). In the latter, the h-h and h-p bonds were taken to have 
transfer matrix elements th-h and th-p equal to 1.2 and 0.9 respectively, so 
that the sum of the transfer integrals remains unchanged after introducing 
alternation. This ratio, which is slightly larger than the experimental value, 
was chosen to magnify the effects of bond alternation and is not expected 
to change the properties qualitatively. We have carried out the simulations 
with At = 0.1, usually for a value of the projection parameter, f3 ^ 2.0. This 
corresponds to a Trotter error of ~ 10^^. We have employed ~ 10"^ MC sweeps. 
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allowing 2000 sweeps for equilibration. We have checked for convergence by 
increasing the value of the projection parameter /?, number of MC sweeps and 
decreasing At and found the above parameters to be satisfactory. 

In table 4 we present the ground state energy, energy/site and energy/bond 
of Ceo with and without bond alternation. Compared to its fragments, Ceo 
has the highest energy/bond at any give value of U/t, indicating that it is 
the most frustrated of the systems studied. The ground state energy increases 
with increasing interaction strength but decreases upon introducing bond al- 
ternation. We compare the ground state energies obtained from the SPQMC 
calculations with those from VMC studies in fig. 4 H, |2^. Since the SPQMC 



calculations are not variationally bounded caution should be exercised in in- 
terpreting the PQMC energies vis a vis the VMC energies. However, test 
calculations on small systems indicate that at intermediate U/t, the absolute 
SPQMC errors in energies are comparable with or smaller than VMC errors. 
We shall presently see that the description of the ground state within the 
SPQMC calculation is better than that afforded by the VMC method. 

Ceo perfers a bond alternating structure in the non-interacting limit since 
even here the h-h bonds have larger bond orders than the h-p bonds |30|, 
Using Coulson's relation this translates into shorter h-h distances compared 



to h-p distances p2|, |33[| . In table 5 we present the bond orders of uniform 
and bond alternating Ceo for different values of the Hubbard parameter. We 
also compare these values with the VMC bond orders in fig 5. The VMC 
bond orders show a stronger dependence on U/t. In table 6 we present the 
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difference in the bond orders of the h-h and the h-p bonds for uniform and bond 
alternant Ceo, within VMC and SPQMC schemes. The VMC calculations seem 
to indicate that interactions supress this difference. However, the SPQMC 
bond orders do not fall so sharply with increasing U /t. In one- dimensional 
Peierls-Hubbard systems, it is known that the difference in bond orders of the 
^double' and ^single' bonds goes through a maximum as a function of ?7/t [^ . 
It is reasonble to expect a similar behaviour in the case of Ceo- However, the 
maximum in this case should shift to larger values of U/t due to increased band 
width of the bands in higher dimensions. The SPQMC bond order difference 
registers a downward trend for U/t between 5.0 and 6.0. 

We examine the spin-spin correlations on fuUerene for different interaction 
strengths (table 5, fig. 6) and for different inter-site separations (fig. 7) for con- 
nectivities shown in the pentagonal projection of Ceo shown in fig. 3. The h-h 
spin-spin correlation is strongly antiferromagnetic while the h-p correlation, 
although afm, is weaker. This trend is already manifest even in corannulene, 
an open system. We observe in fig. 6 the evolution of the nearest neighbour 
and longer range correlations with increasing U/t. As with the bond orders, 
increasing interaction strength enhances the difference between the h-h and 
h-p spin-spin correlations. A Hartree-Fock (HF) calculation on Ceo predicts 
an instabilty in the spin structure of the ground state for U/t ^ 2.6||12|. The 



evolution of the SPQMC nearest neighbor and longer range correlations with 
increasing correlation strength is smooth (fig. 6), showing no indication of any 
instability in the vicinity of U/t ^ 3.0. We also study the charge correlations 
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(table 5 and fig. 8) and find no evidence of any instability in the ground state. 
Thus, the properties of the ground state tend to evolve smoothly from the 
non- interacting to the strongly correlated limit. 

The SPQMC procedure allows a more detailed analysis of the bond orders 
and spin-spin correlations, since all the symmetry related estimates are strictly 
equal. We can construct the molecular structure factors (MSF) of Cqq, similar 
to the structure factors in crystals, by projecting the quantities such as bond 
orders and spin correlations on to different irreducible representations of the 
icosahedral point group. The MSF corresponding to the bond orders for the 
ptft. irreducible representation, pr, is given by 

Pr = ^ E Xr(^) (E 4.cnia + h.c.) ■ (36) 

R 

Similarly the MSF corresponding to the spin correlations Sr is given by 

^r=^Exr(^)«4J (37) 

R 

where Xt{R) is the character of the symmetry operator R in the F*'* irreducible 
representation of the group and s^Sp^. is the spin-spin correlation of the z- 
component of the spins at sites i and Ri. 

In fig. 8a and 8b we show the dependence of pr on the irreducible repre- 
sentation r for uniform and bond alternated fullerenes. Since bond orders are 
empirically related to the bond distances, these quantities give the amplitude 
for the nuclear displacement corresponding to the chosen representation F. We 
find that the amplitudes are very small for the Aig and Aiu representations. 
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The former corresponds to a breathing in and out of the fullerene molecule 
while the latter corresponds to diametrically opposed motions of the nuclei re- 
lated by the inversion symmetry. The maximum amplitude is for the Hg mode 
while the amplitudes for the other modes are almost equal in magnitude. This 
clearly shows that the strongest electron-phonon coupling would be to the Hg 
mode. Varma et a/[^ have indeed employed electron-phonon coupling to this 
mode as a mechanism for superconductivity in fullerides. The effect of increas- 
ing the strength of electron correlations is to reduce the amplitude of the Hg 
mode. Thus, the electron-phonon coupling constant decreases with increase in 
the Hubbard parameter. Bond alternated Ceo also shows features very similar 
to uniform Cgo- 

In fig. 9a and 9b, the dependence of the spin structure factors on the irre- 
ducible representation are shown for uniform and bond alternated fullerenes. 
In both the cases, we find large amplitudes for the T2g and Gu projections. 
Thus, the spin structure in both the cases are similar and correspond to a mix 
of these two representations. A geometrical picture of the spin structure with 
a given V can be obtained by placing individual spins, at each site, directed 
along the site displacement vector for vibrations in that mode[^. It is in- 
teresting to see that the spin MSF has an altogether different structure than 
the bond MSF. Thus, for instance, the nuclear geometry resulting from ap- 
plication of pressure will be qualitatively different from the spin orientations 
resulting from the application of a magnetic field. The effect of increase in 
the strength of electron correlations on the spin structure is complex. While 
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the amplitudes of the dominant representations increases with increase in U/t, 
that of a few others decreases and this could be due to an underlying sum rule 
for the chosen total spin, which in the ground state is a singlet. 

5 Summary 

We have studied the Ceo molecule and its molecular fragments using a new, 
symmetrized YersioTi of the PQMC method developed by us. In this method all 
the symmetry related Ising configurations are simultaneously sampled which 
leads to property estimates that reflect the full symmetry of the system. This 
method also has the potential to project on to the lowest state in subspaces of 
different symmetries. We have tested the SPQMC method on small systems. 
We find a clear improvement in the correlation functions when compared to 
the nonsymmetrized estimates. 

We have applied the SPQMC method to molecular fragments of Ceo- We 
find that corannulene, which has a pentagonal unit completely surrounded by 
ethylenic units is the fragment that resembles Ceo most closely. The interior 
sites of corannulenene resemble radialene which is the basic structural motif 
in Ceo- We have obtained the SPQMC ground states of Ceo with and without 
bond alternation. These have been compared with Hartree-Fock and VMC 
results. SPQMC bond orders indicate that interactions tend to enhance bond 
alternation as in the case of one-dimensional Peierls-Hubbard systems. We find 
no evidence of an instability in the spin structure of Ceo with increasing corre- 
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lation strength. We have analyzed the bond orders and spin-spin correlation 
data in terms of their projections on the space of irreducible representations of 
the icosahedral point group. These molecular structure factors are analogous 
to the structure factors of translationally invariant systems and give the am- 
plitudes for distortions which transform as the irreducible representations of 
the point group. The amplitude for the Hg distortion is largest and its mag- 
nitude decreases with increasing U/t. Therefore, it is expected that the Hg 
vibrational mode couples strongly with the electrons in the system for small 
values oiU/t. The spin structure has a bimodal distribution with large weights 
in the and representations, whose amplitudes increase with increasing 
U/t. 

Acknowledgements: We thank Y. Anusooya for help with the exact com- 
putations on small systems and Dr. Biswadeb Dutta for help with the com- 
putational facihty at JNCASR. 
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Table 1. Ground state energy of benzene (in units of t) as a function of U/t 

using exact diagonalization {E^xact), GWF as the trial wavefunction in a 
Gutzwiller variational calculation {Ecutz.), VMC and SPQMC calculations. 
Quantities in parantheses are absolute values of the percentage errors. 



U/t 


Eexact 


Ecutz 


EvMC 


EsPQMC 


0.5 


-7.274 


-7.274 (i0.007) 


-7.274 (i0.007) 


-7.278 (0.05) 


1.0 


-6.601 


-6.596 (0.08) 


-6.599 (0.03) 


-6.611 (0.15) 


2.0 


-5.410 


-5.386 (0.44) 


-5.388 (0.41) 


-5.444 (0.63) 


4.0 


-3.669 


-3.537 (3.60) 


-3.538 (3.57) 


-3.783 (3.11) 



Table 2(a). Spin-Spin correlations (4 < sfsj >) of benzene for U/t — 1.0 and 
4.0, from exact and MC calculations. Quantities in parantheses are absolute 

values of the percentage errors. 







U/t 


'■J 


1 < sjs) > 






Exact Gutz VMC SPQMC 


1.0 
1.0 
1.0 

1.0 


1,1 
1,2 
1,3 
1,4 


0.567 0.565 (0.35) 0.557 (1.76) 0.562 (0.88) 
-0.267 -0.265 (0.75) -0.245 (8.24) -0.264 (1.12) 
0.022 0.016 (27.3) 0.011 (50.0) 0.020 (9.10) 
-0.077 -0.068 (11.7) -0.062 (19.5) -0.074 (3.90) 


4.0 
4.0 
4.0 
4.0 


1,1 
1,2 

1,3 

1,4 


0.778 0.751 (3.47) 0.750 (3.60) 0.751 (3.47) 
-0.435 -0.404 (7.13) -0.395 (9.20) -0.412 (5.29) 
0.140 0.095 (32.1) 0.077 (45.0) 0.120 (14.0) 
-0.188 -0.131 (30.3) -0.113 (39.9) -0.167 (11.2) 
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Table 2(b). Charge-Charge correlations (< n^nj >) of benzene for U/t — 1.0 
and 4.0, from exact and MC calculations. Quantities in parantheses are the 
absolute values of the percentage errors. 







U/t 




< riiTij > 






Exact Gutz VMC SPQMC 


1.0 

1.0 
1.0 
1.0 


1,1 
1,2 
1.3 
1,4 


1.432 1.435 (0.21) 1.443 (0.77) 1.436 (0.28) 
0.816 0.817 (0.12) 0.812 (0.49) 0.813 (0.37) 
0.989 0.989 (i0.05) 0.976 (1.31) 0.989 (i0.05) 
0.960 0.953 (0.73) 0.952 (0.83) 0.958 (0.21) 


4.0 
4.0 
4.0 
4.0 


1,1 
1,2 
1,3 
1,4 


1.222 1.249 (2.21) 1.250 (2.29) 1.250 (2.29) 
0.905 0.911 (0.66) 0.911 (0.66) 0.896 (0.99) 
0.990 0.979 (1.11) 0.979 (1.11) 0.987 (0.30) 
0.989 0.971 (1.82) 0.971 (1.82) 0.985 (0.40) 



Table 3. Ground state energies of pyracylene, fluoranthene and corannulene 





(in units of t) from VMC and PQMC calculations. 


U/t 


Pyracylene 


Fluoranthene 


Corannulene 




EvMC EpQMC 


EvMC EpQMC 


EvMC EpQMC 


0.5 
1.0 
2.0 
4.0 


-17.713 -17.725 
-16.116* -16.158 
-13.235 -13.406 
-8.666 -9.434 


-20.558 -20.565 
-18.730 -18.771 
-15.422 -15.581 
-10.164 -10.892 


-26.3058 -26317 
-24.007 -24.062 
-19.8474 -20.037 
-13.158 -14.063 



* exact energy obtained from full CI calculation in the singlet subspace of 
dimension 2760615 is -16.13 
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Table 4. Ground state energy of uniform and bond alternant Ceo for various 
values of U from SPQMC calculations. 



u/t 


uniform 


bond alternant 




E 


E/site 


E/bond 


E 


E/site 


E/bond 


1.0 


-79.0307 


-1.317 


-0.878 


-81.5973 


-1.360 


-0.907 


2.0 


-66.6284 


-1.110 


-0.740 


-69.2986 


-1.155 


-0.770 


3.0 


-56.0459 


-0.934 


-0.623 


-58.8326 


-0.981 


-0.654 


4.0 


-47.3407 


-0.789 


-0.526 


-50.2902 


-0.836 


-0.559 


5.0 


-40.4724 


-0.675 


-0.450 


-43.5544 


-0.726 


-0.484 


6.0 


-35.6031 


-0.593 


-0.396 


-38.4948 


-0.642 


-0.428 



Table 5. On-site and nearest neighbour SPQMC charge-charge correlations 
(< niUj >) and spin-spin correlations (4 < >) and bond orders of 
uniform and bond alternant Cqq for different values of U/t. 







u 




uniform 


bond alternant 






h-h 


h-p 


diag 


h-h 


h-p 


diag 


1.0 


bo 


0.598 


0.472 




0.708 


0.408 






s-s 


-0.212 


-0.130 


0.549 


-0.294 


-0.097 


0.550 




c-c 


0.816 


0.902 


1.150 


0.783 


0.926 


1.149 


2.0 


bo 


0.589 


0.461 




0.700 


0.396 






s-s 


-0.249 


-0.148 


0.559 


-0.345 


-0.109 


0.603 




c-c 


0.869 


0.915 


1.401 


0.814 


0.937 


1.397 


3.0 


bo 


0.575 


0.442 




0.687 


0.374 






s-s 


-0.294 


-0.168 


0.650 


-0.406 


-0.121 


0.655 




c-c 


0.888 


0.928 


1.350 


0.840 


0.947 


1.345 


4.0 


bo 


0.554 


0.415 




0.667 


0.345 






s-s 


-0.350 


-0.192 


0.702 


-0.473 


-0.135 


0.708 




c-c 


0.905 


0.939 


1.298 


0.863 


0.912 


1.292 



29 



Table 6. Difference in the bond orders of the h-h and h-p bonds of Ceo from 



VMC and SPQMC calculations. 





uniform 


bond alternant 




VMC SPQMC 


VMC SPQMC 


1.0 
2.0 
3.0 
4.0 


0.126 0.126 
0.121 0.128 
0.121 0.133 
0.112 0.139 


0.482 0.300 
0.486 0.304 
0.436 0.313 
0.379 0.322 
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Figure Captions 

Figure 1. Sample Ising configurations (a) and (b) for the 4-site Hubbard 
ring witli 4 time slices which are related by C4 symmetry. 

Figure 2. Structures of (a) Pyracylene, (b) Fluoranthene and (c) Corannu- 
lene. 

Figure 3. (a) Structure of Radialene and (b) Pentagonal projection of Ceo- 

Figure 4. Energy/site of (a) uniform and (b) bond alternated Cgo for various 
values oiU/t. Open circles represent SPQMC data while crosses are from 
the VMC calculation. 

Figure 5. Bond orders for the h-h and h-p bonds from SPQMC calculations 
(crosses and open triangles) and VMC calculations (open squares and 
open circles) for (a) uniform and (b) bond alternated Cqq for various 
values oiU/t. 
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Figure 6. Inter-site spin-spin correlations of (a) uniform and (b) bond alter- 
nant Ceo ioY 1-2 (filled squares), 1-11 (open circles), 1-3 (asterisks), 1-6 
(crosses) and 1-7 (open triangles) site connectivities, for different values 
of U/t. 

Figure 7. Inter-site charge-charge correlations of (a) uniform and (b) bond 
alternant Ceo for 1-2 (open squares), 1-11 (open circles), 1-3 (asterisks), 
1-6 (crosses) and 1-7 (open triangles) site connectivities, for different 
values oiU/t. 

Figure 8. pr for (a) uniform and (b) bond alternant Ceo for the irreducible 
representations of the icosahedral point group, ior U/t — 1.0 (asterisks), 
U/t = 2.0 (open circles), U/t = 3.0 (triangles), U/t = 4.0 (crosses), 
U/t — 5.0 (inverted triangles) and U/t — 6.0 (open squares). 

Figure 9. Sr for (a) uniform and (b) bond alternant Ceo for the irreducible 
representations of the icosahedral point group, ior U/t — 1.0 (asterisks), 
U/t = 2.0 (open circles), U/t — 3.0 (triangles), U/t = 4.0 (crosses), 
U/t — 5.0 (inverted triangles) and U/t — 6.0 (open squares). 
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